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Abstract — This paper presents a new hybrid particle swarm 
optimization (HPSO) method for solving multi-objective real 
power optimization problem. The objectives of the 
optimization problem are to minimize the losses and to 
maximize the voltage stability margin. The proposed method 
expands the original GA and PSO to tackle the mixed -integer 
non- linear optimization problem and achieves the voltage 
stability enhancement with continuous and discrete control 
variables such as generator terminal voltages, tap position of 
transformers and reactive power sources. A comparison is made 
with conventional, GA and PSO methods for the real power 
losses and this method is found to be effective than other 
methods. It is evaluated on the IEEE 30 and 57 bus test system, 
and the simulation results show the effectiveness of this 
approach for improving voltage stability of the system. 

Keywords: Hybrid Particle Swarm Optimization (HPSO), real 
power loss, reactive power dispatch (RPD), Voltage stability 
constrained reactive power dispatch (VSCRPD). 

I. Introduction 

Optimal reactive power dispatch problem is one of the 
difficult optimization problems in power systems. The sources 
of the reactive power are the generators, synchronous 
condensers, capacitors, static compensators and tap 
changing transformers. The problem that has to be solved in 
a reactive power optimization is to determine the optimal 
values of generator bus voltage magnitudes, transformer tap 
setting and the output of reactive power sources so as to 
minimize the transmission loss. In recent years, the problem 
of voltage stability and voltage collapse has become a major 
concern in power system planning and operation. To enhance 
the voltage stability, voltage magnitudes alone will not be a 
reliable indicator of how far an operating point is from the 
collapse point [1]. The reactive power support and voltage 
problems are intrinsically related. Hence, this paper formulates 
the reactive power dispatch as a multi -objective optimization 
problem with loss minimization and maximization of static 
voltage stability margin (SVSM) as the objectives. Voltage 
stability evaluation using modal analysis [2] is used as the 
indicator of voltage stability enhancement. The modal 
analysis technique provides voltage stability critical areas 
and gives information about the best corrective/preventive 
actions for improving system stability margins. It is done by 
evaluating the Jacobian matrix, the critical eigen values/vector 
[3] .The least singular value of converged power flow jacobian 
is used an objective for the voltage stability enhancement. 
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It is a non- linear optimization problem and various 
mathematical techniques have been adopted to solve this 
optimal reactive power dispatch problem. These include the 
gradient method [4, 5], Newton method [6] and linear 
programming [7]. The gradient and Newton methods suffer 
from the difficulty in handling inequality constraints. To apply 
linear programming, the input- output function is to be 
expressed as a set of linear functions, which may lead to loss 
of accuracy. Recently, global optimization techniques such 
as genetic algorithms have been proposed to solve the 
reactive power optimization problem [8-15]. Genetic algorithm 
is a stochastic search technique based on the mechanics of 
natural selection [ 1 6] .In GA -based RPD problem it starts with 
the randomly generated population of points, improves the 
fitness as generation proceeds through the application of 
the three operators-selection, crossover and mutation. But 
in the recent research some deficiencies are identified in the 
GA performance. This degradation in efficiency is apparent 
in applications with highly epistatic objective functions i.e. 
where the parameters being optimized are highly correlated. 
In addition, the premature convergence of GA degrades its 
performance and reduces its search capability. In addition to 
this, these algorithms are found to take more time to reach 
the optimal solution. Particle swarm optimization (PSO) is 
one of the stochastic search techniques developed by 
Kennedy and Eberhart [17]. This technique can generate high 
quality solutions within shorter calculation time and stable 
convergence characteristics than other stochastic methods. 
But the main problem of PSO is poor local searching ability 
and cannot effectively solve the complex non-linear equations 
needed to be accurate. Several methods to improve the 
performance of PSO algorithm have been proposed and some 
of them have been applied to the reactive power and voltage 
control problem in recent years [18-20]. Here a few 
modifications are made in the original PSO by including the 
mutation operator from the real coded GA. Thus the proposed 
algorithm identifies the optimal values of generation bus 
voltage magnitudes, transformer tap setting and the output 
of the reactive power sources so as to minimize the 
transmission loss and to improve the voltage stability. The 
effectiveness of the proposed approach is demonstrated 
through IEEE-30and IEEE-57 bus system. 
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II PROBLEM FORMULATION 

Power systems are expected to operate economically 
(minimize losses) and technically (good stability). Therefore 
reactive power optimization is formulated as a multi-objective 
search which includes the technical and economic functions. 

A. Economic function: 

The economic function is concerned mainly to minimize 
the active power transmission loss and it is stated as, since 
reduction in losses reduces the cost. 

MinP =/(*-*.)= 2>^ 2 +V/-2V,V,cos^) 
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where f(x l ,x 2 ) denotes the active power loss function of 

the system. 

V G is the generator voltage (continuous) 

T k is the transformer tap setting (integer) 

Q c is the shunt capacitor/ inductor (integer) 

V L is the load bus voltage 

Q G is the generator reactive power 

k = (i,j)J e N B ,J e N itgi is the conductance of branch k. 

dy is the voltage angle difference between bus I &j 

P Gi is the injected active power at bus i 

P Di is the demanded active power at bus i 

Gy is the transfer conductance between bus i and j 

By is the transfer susceptance between bus i and j 
Q Gi is the injected reactive power at bus i 
Q Di is the demanded reactive power at bus i 
N e is the set of numbers of network branches 
N P q is the set of number of PQ buses 
N b is the set of numbers of total buses 

Nj is the set of numbers of buses adjacent to bus i 
(including bus i ) 
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N is set of numbers of total buses excluding slack bus 

N c is the set of numbers of possible reactive power 

source installation buses 
N, is the set of numbers of transformer branches 

Si is the power flow in branch / the subscripts 
'min' and "max" in Eq. (2-7) denote the corresponding 
lower and upper limits respectively. 

B. Technical function: 

The technical function is to minimize the bus voltage 
deviation from the ideal voltage and to improve the voltage 
stability margin (VSM) and it is stated as 

Max (VSM=max (minleig (jacobi)) (8) 

where jacobi is the load flow jacobian matrix , eig (jacobi) 
returns all the eigen values of the Jacobian matrix, 
min(eig( Jacobi)) is the minimum value of eig (Jacobi) , max 
( min ( eig (Jacobi))) is to maximize the minimal eigen value in 
the Jacobian matrix. 

IE. PARTICLE SWARM OPTIMIZATION (PSO) 

A. Overview: 

PSO is a population based stochastic optimization 
technique developed by Kennedy and Eberhart [17]. A 
population of particles exists in the n-Dimensional search 
space. Each particle has a certain amount of knowledge, and 
will move about the search space based on this knowledge. 
The particle has some inertia attributed to it and so it will 
continue to have a component of motion in the direction it is 
moving. It knows where in the search space, it will encounter 
with the best solution. The particle will then modify its 
direction such that it has additional components towards its 
own best position, pbest and towards the overall best 
position, gbest. The particle updates its velocity and position 
with the following Equations (9) to (1 1) 

Vf =W*V? +q *ran® x *ipbe^t-Sf)+Q *ran4) 2 *(gbesjt-S*)(9) 
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: Velocity of particle / at the iteration k + 1 
: Velocity of particle / at the iteration k 
: Position of particle / at the iteration k + 1 
: Position of particle / at the iteration k 

: Constant weighting factor related to pbest 
: Constant weighting factor related to gbest 



rand( ) , : Random number between and 1 
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rand ( ), : Random number between and 1 
pbest i : pbest position of particle i 

gbest : : gbest position of swarm 
Usually the constant weighting factor or the acceleration 
coefficients C l ,C 2 = 2 , control how far a particle moves in 
a single iteration. The inertia weight' W is used to control 
the convergence behavior of PSO. Suitable selection of the 
inertia weight provides a balance between global and local 
exploration and exploitation of results in lesser number of 
iterations on an average to find a sufficient optimal 
solution. In the PSO method, there is only one population 
in an iteration that moves towards the global optimal point. 
This makes PSO computationally faster and the 
convergence abilities of this method are better than the 
other evolutionary computation techniques such as GA. 

B. Proposed Algorithm: 

The main drawback of the PSO is the premature 
convergence. During the searching process, most particles 
contract quickly to a certain specific position. If it is a local 
optimum, then it is not easy for the particles to escape from it. 
In addition, the performance of basic PSO is greatly affected 
by the initial population of the particles, if the initial population 
is far away from the real optimum solution . A natural evolution 
of the PSO can be achieved by incorporating methods that 
have already been tested in other evolutionary computation 
techniques. Many researchers have considered incorporating 
selection, mutation and crossover as well as differential 
evolution into the PSO algorithm. The main goal is to increase 
the diversity of the population by: preventing the particles 
to move too close to each other and collide, to self-adapt 
parameters such as constriction factor, acceleration constants 
or inertia weight. As a result, hybrid versions of PSO have 
been created and tested in different applications. In the 
proposed approach, mutation which is followed in genetic 
algorithm is carried out. Mutation is one of the effective 
measures to prevent loss of diversity in a population of 
solution, which can cover a greater region of the search 
space. Hence in this algorithm the addition of mutation into 
PSO will expand its global search space, add variability into 
the population and prevent stagnation of the search in local 
optima. 

The mutation operator works by changing a particle 
position dimension 



using' niutate{Sj ) - 



S, + delta(iter, U -S,):rb = l 
S, + delta(iter, S, ■ - L) : rb = 



(12) 



Where iter is the current iteration number, 
U is the upper limit of variable spac 
L is the lower limit of variable space 
rb is the randomly generated bit 
delta (iter, y) return a value in the range [0: y] 
It provides a balance between adding variability and allowing 
the particles to converge. Hence in this method it reduces 
the probability of getting trapped into local optima. 
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C. HPSO Algorithm Procedure: 

Step 1 : Initialization of the parameters 

Step 2: Randomly set the velocity and position 
of all the particles. 

Step 3: Evaluate the fitness of the initial 
particles by conducting Newton-Raphson power flow 
analysis results. pbest of e ach particle is set to initial 
position. The initial best evaluation value among the 
particles is set to gbest. 

Step 4: Change the velocity and position of the particle 
according to the equations (9) to (1 1). 

Step 5: Select the best particles come into mutation 
operation according to (12). 

Step 6: If the position of the particle violates the limit 
of variable, set it to the limit value. 

Step 7: Compute the fitness of new particles. If the 
fitness of each individual is better than the 

previous pbest; the current value is set to 

pbest value. If the best pbest is better than 

gbest, the value is set to be gbest. 

Step 8: The algorithm repeats step 4 to step 7 until the 
convergence criteria is met, usually a sufficiently good 
fitness or a maximum number of iterations. 

IV .HPSO IMPLEMENTATION OF THE OPTIMAL 
REACTIVE POWER DISPATCH PROBLEM: 

When applying HPSO to solve a particular optimization 
problem, two main issues are taken into consideration namely: 

(i) Representation of the decision variables and 

(ii) Formation of the fitness function 
These issues are explained in the subsequent section. 

A. Representation of the decision variables 

While solving an optimization problem using HPSO, each 
individual in the population represents a candidate solution. 
In the reactive power dispatch problem, the elements of the 
solution consists of the control variables namely; Generator 
bus voltage (V .), reactive power generated by the capacitor 
(Q ci ), and transformer tap settings (t k ). Generator bus voltages 
are represented as floating point numbers , whereas the 
transformer tap position and reactive power generation of 
capacitor are represented as integers. With this 
representation the problem will look like the following: 

098897010 1U92696.5100CE321G.15I 

ff J$ ^ A ii %i fit fit fii 

B. Formation of the fitness function 

In the optimal reactive power dispatch problem, the 
objective is to minimize the total real power loss while 
satisfying the constraints (14) to (20). For each individual, 
the equality constraints are satisfied by running Newton- 
Raphson algorithm and the constraints on the state variables 
are taken into consideration by adding penalty function to 
the objective function. With the inclusion of the penalty 
function, the new objective function then becomes, 
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Min F = />,,„ + wEig "» + £ W» + £ 2^, + Z LP ' (13) 

i=i ,=i /=i 

where w, A^ ,K cr K t are the penalty factors for the eigen 

value, load bus voltage limit violation, generater reactive 
power limits violation and line flow limit violation respectively 
.In the above expressions 



solution was obtained with the following parameter setting: 



VP 



K v {v,-vrf if v,>v; 

K v {V t -Vrf if V t <V? 
otherwise 



(14) 



e^ 
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otherwise 



Ik^-^T) 1 ifs t >s? 
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"0 



otherwise 



(16) 



Generally, PSO searches for a solution with maximum fitness 
function value. Hence, the minimization objective function 

given in (17) is transformed into a fitness function (/) to be 
maximized as, 

f = K/F (17) 

where K is a large constant. This is used to amplify (1/F), the 
value of which is usually small, so that the fitness value of 
the chromosomes will be in a wider range. 

V.SIMULATION RESULTS 

In order to demonstrate the effectiveness and robustness 
of the proposed technique, minimization of real power loss 
under two conditions, without and with voltage stability 
margin (VSM) were considered. The validity of the proposed 
PSO algorithm technique is demonstrated on IEEE- 30and 
IEEE-57 bus system. The IEEE 30-bus system has 6 generator 
buses, 24 load buses and 41 transmission lines of which 
four branches are (6-9), (6-10) , (4-12) and (28-27) - are with 
the tap setting transformers. The IEEE 57-bus system has 7 
generator buses, 50 load buses and 80 transmission lines of 
which 17 branches are with tap setting transformers. The real 
power settings are taken from [1]. The lower voltage 
magnitude limits at all buses are 0.95 p.u. and the upper limits 
are 1.1 for all the PV buses, 0.05 p.u. for the PQ buses and the 
reference bus for IEEE 30-bus system. The PSO -based 
optimal reactive power dispatch algorithm was implemented 
using the MATLAB programmed and was executed on a 
Pentium computer. 

Case A : RPD with loss minimization objective 

Here the PSO-based algorithm was applied to identify the 
optimal control variables of the system .It was run with 
different control parameter settings and the minimization 
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:0.4 


c, 




:2 


c 2 




:2 


Maximum generations: 50 


Mutate rate 


:0.1 





Figure 1 illustrates the relationship between the best fitness 
values against the number of generations. 




Figure . 1. Convergence characteristics 

From the figure it can be seen that the proposed algo- 
rithm converges rapidly towards the optimal solution. The 
optimal values of the control variables along with the mini- 
mum loss obtained are given in Table I for IEEE-30 bus sys- 
tem. Corresponding to this control variable setting, it was 
found that there are no limit violations in any of the state 
variables. To show the performance of the HPSO in solving 
this integer nonlinear optimization problem, it is compared to 
the well known conventional, GA &PSO techniques. But in 
HPSO the best solution is achieved. This shows HPSO is 
capable of reaching better solutions and is superior compared 
to other methods. This means less execution time and less 
memory requirements. 

TABLE I 
RESULTS OF PSO-RPD OPTIMAL CONTROL VARIABLES 



Control Variables 


RPD 


VSCRPD 


Vl 


1.0500 


1.054B 


V 2 


1.0504 


1.070 


Vs 


1.0507 


1.004B 


Va 


1.0443 


1.0095 


V tl 


1.1000 


0.9B33 


V 13 


1.0134 


1.05 


Tll 


1.0723 


1.0429 


Tn 


1.1000 


0.92S6 


Tu 


0.942B 


1.0429 


T 3 * 


0.9140 


0.9 


Qcio 


2.5104 


5 


QcLa 


4.9977 





QcLi 


5 





CkL- 


5 


2.B571 


QC3« 


5 


3.5714 


Qcai 


5 


4.2B57 


QC23 


1.5100 


0.7143 


Qcat 


2.5366 


0.7143 


QtL!3 


2.6442 


5 


Loss 


4.5756 


S.2BB7 


VSM 


0.2403 


0.2437 



19 



■A: ACEEE 



ACEEE International Journal on Electrical and Power Engineering, Vol. 1, No. 2, July 2010 



Case B: Multi-objective RPD (RPD including voltage stability 
constraint) 

In this case, the RPD problem was handled as a multi- 
objective optimization problem where both power loss and 
maximum voltage stability margin of the system were 
optimized simultaneously. The optimal control variable 
settings in this case are given in the last column of Table I. To 
maximize the stability margin the minimum eigen value should 
be increased. Here the VSM has increased to 0.2437 from 
0.2403, an improvement in the system voltage stability. For 
IEEE-57 bus system the minimum power loss obtained is 
25.6665 MW.The VSM has increased to 0.1568 from 0. 1456. 
To determine the voltage security of the system, contingency 
analysis was conducted using the control variable setting 
obtained in case A and case B. The eigen values 
corresponding to the four critical contingencies are given in 
Table II. From this result it is observed that the eigen values 
has increased appreciably for all contingencies in the second 
case. This improvement in voltage stability was achieved 
because of the additional objective included in the RPD 
problem in the base case condition. This shows that the 
proposed algorithm has helped to improve the voltage stability 
of the system. To analyze the simulation results it has been 
compared with other optimization methods. Table III 
summarizes the minimum power loss obtained by these 
methods for the IEEE-30 bus system. 



VSM 



TABLE II 
UNDER CONTINGENCY STATE 



Nd 


ContinEmcy 


V5M 


1 

2 

3 
4 


28-27 
4-12 
9-10 

6-9 


RPD 


V5CRPD 


0.13*7 
0.1660 
0.1045 
0.1023 


0.135S 
0.169S 
0.1091 

0.1055 



TABLE III 

COMPARISON OF OPTIMAL RESULT OBTAINED BY 

DIFFERENT METHODS FOR IEEE-30 BUS SYSTEM 



Nd 


Method 


Minimum 

power 
lo5s{MW) 


I 
2 
3 

4 


F v nlnti nTTflfv jifin pramTniTiRf 1 41 

Gradient based optimization^! 4) 
Genetic AlEorithmO. 3) 
Particle swarm optimization^! S) 


4.5613 

5.7360 

5.32 

4.61 



CONCLUSION 

This paper presents a hybrid particle swarm optimization 
algorithm approach to obtain the optimum values of the 
reactive power variables including the voltage stability 
constraint. The effectiveness of the proposed method for 
RPD is demonstrated on IEEE-30 and IEEE-57 bus system 
with promising results. Simulation results show that the HPSO 
based reactive power optimization is always better than those 
obtained using conventional, GA and simple PSO methods. 
From this multi-objective reactive power dispatch solution 
the application of HPSO leads to global search with fast 
convergence rate and a feature of robust computation. Hence 
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from the simulation work, it is concluded that PSO performs 
better results than the conventional methods. 

REFERENCES 

[I] C.A. Canizares, A.C.Z.de Souza and V.H. Quintana, 
"Comparison of performance indices for detection of proximity to 
voltage collapse," vol. 11. no. 3 , pp. 1441-1450, Aug 1996. 

[2] B.Gao ,G.K Morison P.Kundur , 'voltage stability evaluation 
using modal analysis ' Transactions on Power Systems ,Vol 7, No 
.4 .November 1992 [9]. 

[3] Taciana .V. Menezes, Luiz .C.P.da silva, and Vivaldo F.da Costa," 
Dynamic VAR sources scheduling for improving voltage stability 
margin," IEEE Transactions on power systems, vol 18,no.2 ,May 

2003 

[3] O.Alsac, and B. Scott, "Optimal load flow with steady state 

security", IEEE Transaction. PAS -1973, pp. 745-751. 

[4] Lee K Y ,Paru Y M , Oritz J L -A united approach to optimal 

real and reactive power dispatch , IEEE Transactions on power 

Apparatus and systems 1985: PAS-104 : 1147-1153 

[5] A.Monticelli , M .V.F Pereira ,and S. Granville , "Security 

constrained optimal power flow with post contingency corrective 

rescheduling" , IEEE Transactions on Power Systems :PWRS-2, 

No. 1, pp. 175-182. ,1987. 

[6] Deeb N, Shahidehpur S.M, Linear reactive power optimization 

in a large power network using the decomposition approach. IEEE 

Transactions on power system 1990: 5(2) : 428-435 

[7] D. Devaraj, and B. Yeganarayana, "Genetic algorithm based 

optimal power flow for security enhancement", IEE proc- 

Generation. Transmission and. Distribution; 152, 6 November 2005. 

[8]- Deb, K. (201): Multi - objective optimization using 

evolutionary algorithms 1st ed. (John Wiley & Sons, Ltd.). 

[9] Q. H. Wu and J. T. Ma, "Power System Optimal Reactive 

Power Dispatch Using Evolutionary Programming", IEEE Trans. 

on Power Systems, Vol. 10, No. 3, pp. 1243-1249, August 1995. 

[10] K.Iba,' Reactive power optimization by genetic algorithm ," 

IEEE Trans.power syst.vol.9,pp.685-692,May 1992 . 

[I I] D.Devaraj, "Improved genetic algorithm for multi - objective 
reactive power dispatch problem," European Transactions on 
Electrical Power, 2007; 17; 569- 581. 

[12] A.J.Urdaneta, J.EGomez, E.Sorrentino, L.Flores, and R.Diaz, 

'A hybrid genetic algorithm for optimal reactive power planning 

based upon successive linear programming," IEEE Trans.power 

syst, vol.14, pp. 1292-1298, Nov. 1999 

[13] Y.Liu,L.Ma, and Zhang ," GA/SA/TS Hybrid algorithm for 

reactive power optimization, ' in Proc. IEEE Power Eng.Soc. Summer 

Meeting,vol.l, pp. 245- 249,July 2000. 

[14] LL.LaiJ.T.Ma /'Application of evolutionary programming 

to reactive power planning -Comparison with non linear 

programming approach," IEEE Trans.Power syst.,vol,12,pp.l98- 

206,Feb.l997 

[15] Kwang Y. Lee and Frank F.Yang, "Optimal Reactive Power 

Planning Using evolutionary Algorithms: A Comparative study for 

Evolutionary Strategy, Genetic Algorithm and Linear Programming", 

IEEE Trans, on Power Systems, 

Vol. 13, No. 1, pp. 101- 108, February 1998. 

[16] D Goldberg, "Genetic algorithms in search, optimization and 

machine learning", Addison-Wesley, 1989. 

[17] J. Kennedy and R. Eberhart, "Particle swarm optimization", 

Proceedings of the IEEE International Conference on Neural 

Networks, Vol. IV, pp. 1942- 1948, 

1995. 



20 



^VACEEE 



ACEEE International Journal on Electrical and Power Engineering, Vol. 1, No. 2, July 2010 

[20]H.Yoshida,K.Kawata,Y.Fukuyama,"A Particle Swarm 

[18 I S.Durairai, P.S.Kannan and D.Devarai , International journal .. . .. ~ .. , ,. . , ... 

J J J optimization tor reactive power and voltage control considering 

of emerging electric power systems, The Berkeley Electronic Press ,. .. ^ „ __„„ 

00 . . , voltage security assessment, IEEE Transactions on power 

.vol 4, issue 1 ,2005,article 1082 pp. 1-15. . , . - ... , ~ nnn 

' „, ' ,„ ' . . T , ,,, „-, ., . systems,vol.l5.no.4,November2000. 
[19J I.G.Vlachogiannis , K.Y.Lee , Contribution ot generation to 

transmission system using parallel vector particle swarm 

optimization", IEEE Transactions on power 

systems,20(4),2005, 1765-1774. 



21 
© 2010 ACEEE JL* A p C C C 

DOI:01.ijepe.01.02.04 AT/AOCCC 



